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We investigate, via Monte Carlo simulations, the phase structure of a system of closed, non- 
intersecting but otherwise non-interacting, loops in 3 Euclidean dimensions. The loops correspond 
to closed trajectories of massive particles and we find a phase transition as a function of their 
mass. We identify the order parameter as the average length of the loops at equilibrium. This 
order parameter exhibits a sharp increase as the mass is decreased through a critical value, the 
behaviour seems to be a cross-over transition. We believe that the model represents an effective 
description of the broken-symmetry sector of the 2+1 dimensional abelian Higgs model, in the 
extreme strong coupling limit. The massive gauge bosons and the neutral scalars are decoupled, 
and the relevant low-lying excitations correspond to vortices and anti-vortices. The functional 
integral can be approximated by a sum over simple, closed vortex loop configurations. We present 
a novel fashion to generate non-intersecting closed loops, starting from a tetrahedral tessellation 
of three space. The two phases that we find admit the following interpretation: the usual Higgs 
phase and a novel phase which is heralded by the appearance of effectively infinitely long loops. We 
compute the expectation value of the Wilson loop operator and that of the Polyakov loop operator. 
The Wilson loop exhibits perimeter law behaviour in both phases implying that the transition 
corresponds neither to the restoration of symmetry nor to confinement. The effective interaction 
between external charges is screened in both phases, however there is a dramatic increase in the 
polarization cloud in the novel phase as shown by the energy shift introduced by the Wilson loop. 
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I. INTRODUCTION 

Although it is one of the simplest gauge theories, the 
abelian Higgs model is of substantial theoretical interest 
[TJ[6]. It corresponds to scalar electrodynamics consisting 
of a charged scalar field and a neutral vector field which 
is the gauge boson of a [7(1) local gauge symmetry. It 
can be defined in any number of space-time dimensions. 

In 1+1 dimensions, the general absence of spontaneous 
symmetry breaking (2] poses a puzzle as to the manifesta- 
tion of the U(l) symmetry for the naively spontaneously 
broken sector. Indeed, topological vortices play the role 
of instantons and give rise to tunnelling transitions which 
end up disordering the vacuum [3] . The symmetry is ac- 
tually restored; however, a U(l) gauge theory in 1+1 di- 
mensions is classically linearly confining. Consequently, 
charged states are hidden in neutral bound states. 

In 2+1 dimensions, the compact version of the theory 
behaves quite differently than the non-compact version. 
A U{1) gauge theory can be thought of as a theory either 
with gauge group U(l) living on the compact manifold 
S 1 , or with gauge group R 1 , (the real numbers under ad- 
dition) living on the non-compact manifold R 1 . The com- 
pact theory, in the unbroken phase, shows linear confine- 
ment of charges, instead of the classically expected log- 
arithmic potential [TJ, due to magnetic monopoles which 
act as instantons. The actual details of the mechanism 
of this confinement are rather complicated and we will 
not describe them here. In the non-compact case, mag- 
netic monopoles do not exist; hence the expression of the 



symmetry should be along more traditional lines: either 
the symmetry is manifest with a logarithmic potential 
between charged particles, or it is spontaneously broken 
and the interaction is screened. In principle, the theory 
could even be linearly confining for fractionally charged 
external sources. 

The classically spontaneously broken sector of the non- 
compact theory in 2+1 dimensions will be of interest in 
this article. Here, the theory has topological solitons, 
Nielsen-Olesen [8] vortices, which tend to disorder the 
vacuum. Vortex lines in a type II superconductor [5] 
are examples of physical phenomena which are well de- 
scribed by such solitons. Looking at the 3-dimensional 
Euclidean version of the theory, these vortices extend to 
tubes of quantized magnetic flux. For these configura- 
tions to have relevance to the functional integral, finitc- 
ness of the action requires that they form closed loops. 
The contribution of such closed vortex loops to the ex- 
pectation value of the Wilson loop[5] was computed, at 
strong coupling, in a heuristic semi-classical analysis by 
Samuel [TU]. There it was proposed that, if the vortices 
are light enough, they should effectively condense, giv- 
ing rise to a novel phase, what was called the "spaghetti 
vacuum." What this means is that contributions to the 
Euclidean functional integral come preponderantly from 
configurations which are full of vortex loops. It was fur- 
ther deduced that there should be a logarithmic potential 
induced between external charges. Such a potential is in 
fact confining: it takes an infinite amount of energy to 
move two particles infinitely far away from each other, al- 
though it is not linearly confining. A phase transition go- 
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ing from the standard symmetry broken phase to a novel 
phase corresponding to a disgorging of vortex tubes into 
the vacuum has also been proposed by Einhorn and Savit 
[TT] in their study of the lattice abelian Higgs model. 

In this paper we study, by means of Monte Carlo simu- 
lations on a lattice [12] , a discretized, effective version 
of the abelian Higgs model. This amounts to the study of 
a gas of loops on a lattice with Boltzmann weight corre- 
sponding to the total length of the loops. We find indeed 
that there is a rather sharp transition from small av- 
erage loop length to a configuration with an effectively 
infinite loop, the average loop length showing a remark- 
able increase. This kind of transition is very reminiscent 
of percolation type transitions [13]. In the Abelian Higgs 
model interpretation the transition is from the standard 
Higgs phase to a novel phase characterized by saturation 
of the functional integral by configurations that are filled 
with vortex flux loops. We do not however find the corre- 
sponding classical logarithmic potential induced between 
external charges [TU] . External charges are still screened; 
however, we find that the energy of the screening cloud 
increases dramatically. 

II. EFFECTIVE ABELIAN HIGGS MODEL AT 
STRONG COUPLING 

The abelian Higgs model is described by the La- 
grangian density 

C = -\F^ + \d^{D^)* - ^(|0| 2 - r? 2 ) 2 , (1) 

where is a complex scalar field, is a U(l) gauge field, 
0^=8^- ieA^, = <9 P A, - d v A tl (ju, v = 0, 1, 2) and 
A, e and r\ are taken to be positive constants. This theory 
undergoes spontaneous symmetry breaking appended by 
the Higgs mechanism yielding a perturbative spectrum of 
a massive vector boson with mass M = erj and a neutral 
scalar boson with mass m = \/2\ r\. 

Additionally, the theory contains vortex solitons of 
quantized magnetic flux in this sector. Their mass be- 
haves like /i = rf x /(2A/e 2 ) [H], where f(2X/e 2 ) is a 
function that satisfies /(l) = 1, but can take any posi- 
tive value as a function of A/e 2 . We can take the strong 
coupling limit A, e — > 00 while keeping 77 and A/e 2 fixed. 
This decouples the perturbative excitations, m, M — > 00, 
leaving only the vortices as the effective excitations. As 
was shown in [8], in this limit, the size of the vortices van- 
ishes and their world lines resemble perfect, fundamen- 
tal strings. We will only study the abelian Higgs model 
in the description afforded by this effective model. The 
phase structure of the effective model must be the same 
as that of the original abelian Higgs model sufficiently 
deep in the strong coupling regime. Thus our results will 
shed light on the asymptotic region of the strong coupling 
limit of the abelian Higgs model. 

In the lowest approximation, neglecting gradient ener- 
gies due to curvature, the action is given by // x L for 



a closed loop of length L. The Euclidean vacuum-to- 
vacuum amplitude is obtained by functionally integrat- 
ing over field configurations that correspond to the fol- 
lowing Euclidean time histories: they are the classical 
vacuum configuration at the initial time, contain a num- 
ber of virtual vortex anti-vortex pairs at intermediate 
steps, and revert back to the classical vacuum configu- 
ration at the final Euclidean time. Thus, in this limit, 
the abelian Higgs model is equivalent to a gas of massive 
particles that carry a conserved flux; these particles are 
non-interacting except when they are in close proximity. 
All other excitations and their interactions are negligible. 
Thus the functional integral is evaluated by integrating 
over field configurations corresponding to closed vortex 
loops [10], indeed a 3 dimensional loop gas. We will calcu- 
late this integral by a numerical Monte Carlo simulation 
on a lattice discretized version of this effective theory. 



III. LATTICE LOOPS 

On the lattice, it is not straightforward how to con- 
struct closed loops. We construct closed loops by start- 
ing with a tessellation of Euclidean 3-space with (non- 
regular) tetrahedra. To generate this tessellation, we 
start with a body-centered cubic (bec) lattice in a box of 
size N = N^N r . Joining the central vertex in each cube 
with its corners fills each cube with 6 identical pyramids 
with square bases given by the faces of each cube. Split- 
ting each pyramid in half yields two irregular tetrahedra 
and the desired tessellation. To define the splitting, we 
start with the cube with one vertex at the origin, extend- 
ing into the positive octant. We cut each face from the 
origin to the opposite diagonal corner in the (x,y) (y,z) 
and the (z,x) plane respectively. Then we translate this 
scheme throughout the lattice. This converts each pyra- 
mid into two (identical) non regular tetrahedra, giving a 
total of 12 tetrahedra in each cube. All points inside the 
box fall into one tetrahedron or another, except for the 
set of measure zero which resides on the surfaces of the 
tetrahedra. Therefore we have filled space with tetrahe- 
dra. 

Loops are generated by distributing the three cube 
roots of unity over the vertices of the tessellation. A 
given triangular face is associated with an oriented length 
of vortex tube piercing it and going to the center if the 
change of phase about the triangle corresponds to ±2tt, 
using the right hand rule. If a triangular face of a given 
tetrahedron has the cube roots of unity distributed on 
the vertices so that one unit of flux enters the tetrahe- 
dron, with a little reflection it is easy to see that as- 
signing any cube root of unity to the fourth vertex of 
the tetrahedron necessarily defines one unit of flux ex- 
iting the tetrahedron through another of its triangular 
faces, passing from the center of the original tetrahedron 
to the center of the neighbouring one. But since space 
is filled with tetrahedra this exiting flux tube enters the 
neighbouring tetrahedron, and repeating the argument it 
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must exit this tetrahedron, entering a third tetrahedron, 
and so on. It is topologically impossible for the vortex 
line to end; it must ultimately close on itself, forming a 
closed loop, since the lattice is made up of only a finite 
number of tetrahedra. The resulting configuration is a 
system of closed vortex loops which are by construction 
non-intersecting, since it is also topologically impossible 
for two vortices to enter the same tetrahedron. A simi- 
lar scheme was originally implemented on a cubic lattice 
[15j for cosmic strings. There, one could have two vortex 
lines entering a cube and two exiting it, although it was 
impossible to resolve the path of the vortex lines inside 
the cube. Our tetrahedral dissection of the cube resolves 
this ambiguity. 

The actual Euclidean geometrical length of the loop 
will depend on the explicit trajectory that the loop takes 
since the distances between the centers of neighbour- 
ing tetrahedra are not all the same. For a long loop, 
these geometrical factors average out, simply giving a 
renormalization of the value of fi, which includes the 
lattice spacing. The corresponding effective action is 
Scfi oc (fi x Lt) where the Lt is the number of trian- 
gles through which the loop passes, where by abuse of 
notation we use the same symbol fi to represent the mass 
times the lattice spacing. We will call Lt the length of 
the loop. The shortest closed loop has Lt = 4 while the 
maximum is LT.max = 12N^N T . 

IV. MONTE CARLO SIMULATIONS 

Our simulations are performed on a bcc cubic lattice 
with N s = 100, N T = 100 and n from to 1.5 using 
Monte Carlo simulations [TT1 [18] . We begin with an ini- 
tial arbitrary configuration of closed loops. Then we use 
the standard Metropolis algorithm to generate an ensem- 
ble of configurations which follow the Boltzmann distri- 
bution with weight given by (e~' tLr ). 
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FIG. 1: (color online) Total length of loops as a func- 
tion of time for 100 3 bcc lattice. From top to bottom, 
fj, = 0,0.1,0.15,0.3,0.9,1.5. Insert displays fi = 0.3, for two 
different initial configurations. 
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A. Thermalization 

In Fig. [T] on a semi logarithmic scale, we show the 
convergence of the total length of loops Lt versus Monte 
Carlo time (updates) for several values of /i and N T — 
100. Our unit of time corresponds to one complete up- 
date of each site of the 100 3 bcc lattice. The equilib- 
rium state does not depend on the initial state, but it is 
strongly dependent on fi. The average total length and 
the absolute fluctuations grow as \x is decreased, but the 
relative fluctuations diminish. 



B. Numerical evidence for a change in phase 

In Fig. [2j we show the expectation value of the total 
length of loops (L T ) as a function of inverse /iona log- 
arithmic scale. We see that there is a dramatic change 



FIG. 2: (color online) Expectation value of the total length of 
loops as a function of fi, < /i < 1.5, with the insert focussing 
on the region of the transition. 

in the curve around fi = 0.15 indicating a transition in 
the system. We define the total density of the loops p 
as the ratio of the computed Lt to £<r imax . The transi- 
tion corresponds to the appearance of effectively infinite 
loops in the simulation. If the simulation could be done 
in infinite space, at the transition a truly infinite loop 
would appear. Infinitely long loops \TE\ in a finite volume 
are operationally defined as those loops having a length 
L, much longer than that they would normally have if 
they corresponded to a closed, self-avoiding random walk. 
The size of a self-avoiding random walk on a simple cubic 
lattice behaves as L~ 3 / 5 , closing the loop adds one con- 
straint which should not greatly modify this scaling law. 
Actually it is found that the exponent should be slightly 




FIG. 4: (color online) Density of finite loops as a function of 
the total density for < [i < 1.5. 




FIG. 3: (color online) Snapshots of the closed loops generated 
in equilibrium phase for 100 3 bcc lattice. Identification of the 
axes is arbitrary. (Top): Only the largest closed loops are 
displayed for /i = 0.152, 1400 < L < 3000. (Bottom): Larger 
loops are formed after the transition (see text) as shown for 
H — 0.148; the closed loop shown has L ~ 155500. 



less than 3/5 on a cubic lattice, but the exact value is not 
analytically known [15]. On our 100 3 bcc lattice, with the 
tetrahedral trajectories, it is not clear how the size of self 
avoiding, closed random walks would exactly scale with 
their length. However, taking the cubic lattice result as 
a guide, simple calculation yields 200 5 / 3 ~ 6840, where 
200 is the number of steps from one side to the other 
of our lattice, hopping from the corner of a cube to the 
center, back to the next corner and so on. Therefore we 
treat any loop of length substantially greater than 10000 
as an infinite loop. 

Fig. [3] shows snapshots of closed vortex loops, with pe- 
riodic boundary conditions, generated in the equilibrium 
phase before the transition, for /i = 0.152 (top) where 
only finite closed loops are present, and after the tran- 
sition, for /i = 0.148 (bottom) where larger closed loops 
are formed. For a better visualization, only some loops 
are presented, as otherwise the picture looks completely 
black, filled with vortex loops. There is some theoretical 
understanding of this phenomenon in thermodynamical 



studies of a network of cosmic strings. The authors of 
[5U] have noted that at formation, the density of states 
of cosmic strings in the early Universe is dominated by 
a large loop containing most of the energy with a ther- 
mal distribution of finite, low mass, strings. This den- 
sity distribution was also described by the authors in 
[15] , The number of microstates available to the sys- 
tem is much greater when reorganized as a large number 
of finite loops augmented by one infinite loop. The den- 
sity at which this happens corresponds to the Hagedorn 
temperature and the transition corresponds to the Hage- 
dorn phase transition [16] . In the cosmological situation, 
as the Universe expands the phase space favours config- 
urations where all the strings are chopped off into the 
smallest possible loops. 

In Fig.|3J we graph the density of finite loops. For small 
values of the total density, there are no infinite loops; 
hence the curve is linear with slope 1. We see a dramatic 
transition around p = 0.207 which corresponds to fi = 
0.152. At the transition there is a sudden reorganization 
of the vortex loops into one infinitely long loop and a 
number of finite loops. Remarkably, the increase in the 
total density /length of loops caused by further decreasing 
\i occurs only by appending to the infinitely long loop, 
the density of finite vortex loops remaining essentially 
constant. 



ORDER PARAMETERS 



We want to analyze the nature of the novel phase and 
to study the system around the transition point. For 
that, we turn to the following observables as order pa- 
rameters: the average length of the loops, the Wilson 
loop operator [H] and the Polyakov [7] loop operator. 



5 



A. Length of the average loop 

The average length of the loops, (L), that make up the 
ensemble of equilibrium configurations shows a remark- 
able transition as a function of fj,. Below we plot (L) as a 
function of In /z. We observe that the transition occurs at 



finite closed vortex loops will give a perimeter behaviour 
for the expectation value of Wilson loop operator, which 
means that the charges are screened. In the novel phase, 
however, the infinitely long vortex loops could give a con- 
tribution that has no relation to the perimeter. 




ln(n) 



FIG. 5: (color online) Average length of the loops as a func- 
tion of ln(/i). 

ln(/z) « —1.9 w ln(0.15), exactly at the point where the 
effectively infinite loop appears. Continuing the figure 
to larger values of ln(/i) is not feasible as the equilib- 
rium configuration is not easily obtained. The Monte 
Carlo process of attaining equilibrium is asymptotically 
slowed. 



B. Wilson loop 

The Wilson loop operator corresponds to inserting into 
the system two static, equal but opposite charges g, sep- 
arating them by a distance L for a duration T with 
(T L), and then annihilating them. The expectation 
value of Wilson loop operator is given by 



W(L 7 T) = (Je^tM^X 



(2) 



where the integration is over the rectangular Wilson loop 
(L x T) contour. For our effective model a dramatic 
simplification occurs, § A^dx^ exactly measures 2tt times 
the linking number v of the Wilson loop with the closed 
vortex loops: 



W(L,T) = (e 



(3) 



For large T, W (L, T) ~ e - A ( L ) T , where A (L) = 
lim T ->-oo -(l/T)\n(W(L,T)), the energy shift, is the in- 
teraction energy of the static qq pair separated by a dis- 
tance L [5]. In the usual Higgs phase, we expect that 
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FIG. 6: (color online) Wilson loop W{L = 20, T = 80) as 
a function of 2-nq/e from to -n for fi = 0.3,0.2,0.152 (up- 
pcr)and /i = 0, 0.13, 0.14, 0.147, 0.152. (lower) 

Figures |6]-|8] shows our results for the numerical calcu- 
lation of the potential between qq. In Fig. [6]we start with 
the calculation of the Wilson loop operator W(L, T) for 
< 2-Kq/e < 2tt for various values of /j, in the Higgs phase 
in the upper graph for [i = .3 to the border of the phase 
transiition at /i = .152 and in the lower graph in the novel 
phase for /i = .152 to fj, = 0. We note that the curve 
moves rapidly from the borderline value at [i = 0.152 to 
that deep in the Higgs phase at /i = 0.3, conversely there 
is very little movement from the transition at [i = .152 
to fJL = 0. 

The energy shift at finite T, A(L,T), is analyzed for 
the value fi = .2 in the Higgs phase, shown in the Fig. [7] 
The energy shift should be a periodic function of the 
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FIG. 7: (color online) dA/dL as a function of 2irq/e for 
T = 80, 90, 98. The data points are our numerical re- 
sults for /i = 0.2. The solid lines are the fit of the form 
C(T) sin 2 ((27r (? /e)/2). 



fit. In Fig. [8j the T dependence of C(T) is displayed; 
it is a linear function of l/T. The extrapolation of the 
curve to T — > oo yields C(oo) — 0, i.e. B = 0. This 
means that there is only the perimeter law behavior for 
the Wilson loop operator in the Higgs phase. 

In the novel phase, for fi < 0.152, the Wilson loop, 
as a function of 2irq/e, does not vary greatly with /i. It 
decreases as x = 2-itq/e approaches approximately 0.98 
where it vanishes. This implies that the interaction en- 
ergy of the external charges is so large that our numerical 
analysis is not able to resolve its value, within the resolu- 
tion permitted by our lattice approximation. It does not 
by any means imply confinement. In the novel phase we 
cannot use the simple function sin 2 ((27rc7/e)/2) to give 
the dependence on q. However we can easily see that 
the Wilson loop is independent of the area. In Fig. [9] we 
plot the Wilson loop as a function of L for a loop of size 
(L, 100 — L) ie. for a fixed value of the perimeter, and 
for a fixed value 2nq/e = .4. It is evident that the value 
of the Wilson loop does not vary. 
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FIG. 8: (color online) C(T) as a function of l/T, fi = .2. 
The dotted line is the linear fit 



external charge q; the expected form is ~ sin 2 ((27rg/e)/2) 
[21) . A perimeter law then would imply 

A(L, T) = A(T) ■ sin 2 ((2^/e)/2)(L + T)/T, (4) 

while an area law would give 

A(L,T) = B(T) ■ sin 2 ((27r(7/e)/2)(L x T)/T. (5) 

In general, we allow a sum of the two behaviours and 
A(T) -> A and B(T) -> B independent of T as T — > oo. 
Then 

dA(L, T)/dL — > (A/T + B)sin 2 ((27rt7/e)/2) (6) 
= C{T) sin 2 ((27r 9 /e)/2) (7) 

should be independent of L. The dots correspond to our 
numerical simulation; the solid lines correspond to the 
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FIG. 9: (color online) Wilson loop W(L, T = 100 - L) as a 
function of area for fi — .13. 

To compare the novel phase with the Higgs phase, we 
can look at the distribution of the linking number, is, of 
the Wilson loop. The histogram for the distribution of 
the linking number are given in Figures [10] and [TT] for 
fi = .13 and fi = .17 respectively. Such histograms 
were calculated for different values of fi on either side of 
the transition and for different sizes of the Wilson loop. 
In Fig. 12 we plot the variance a 2 , of the histograms 
as a function of the perimeter, L + T for L — 40 and T 
varying. Clearly we find a perimeter law for the variance. 

The expectation value of Wilson loop is constructed 
from the histograms of linking number of the Wilson 
loop. We compute e l ( 27ri 3/e)xi/ j? or eacn V alue of the link- 
ing number v, weigh that phase with the number of con- 
figurations with that value of v, and then sum over all 
linking numbers. This actually corresponds to calculat- 
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FIG. 10: (color online) Histogram of the Wilson loop linking 
number v for W(L, T — 40, 80) for /i = .13 in the novel phase. 
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FIG. 11: (color online) Histogram of the Wilson loop linking 
number for W(L, T — 40, 80) for /i = .17 in the Higgs phase. 



ing the characteristic function of the probability distri- 
bution function for the linking number [22]. Clearly this 
gives a sum of phases which are distributed over the unit 
circle in a manner depending on the exponent. If the 
a x 2irq/e is of the order of o(l), the phases are essentially 
randomly distributed over the unit circle, and the char- 
acteristic function vanishes. This is seen in Fig. [6] (lower 
graph) where the expectation value of the Wilson loop 
crashes to zero for 2nq/e ~ .6 in the novel phase, when 
the variance suddenly becomes large. This behaviour 
is to be contrasted with that in the Higgs phase where 
the expectation value is a smooth sinusoidal function of 
2-Kq/e £ [0, 27r] as seen in Fig. [6] (upper) for > -152. 

The perimeter law for the variance does translate into 
a perimeter law for the Wilson loop. In Fig. [13] and [14] 
for /i = .2 in the Higgs phase and for fj, = .13 in the 
novel phase, we plot the the log of the expectation of the 
Wilson loop as a function of the perimeter. We see a 
perimeter law for the Wilson loop and a linear behaviour 
of the derivative of the energy shift in both phases. We 
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FIG. 12: (color online) Plots of the variance as a function of 
the perimeter for /i = .13, .2, .3 and .6. The dashed lines are 
a linear fit. 



do our analysis for a fixed value of 2nq/e = .4. 
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FIG. 13: (color online) Plot of -ln(W) for fi = .13. The 
dashed lines are a linear fit, we find — ln(W) = .009280(±4 x 
10" 6 ) • (L + T) + .00350(±5 x 10~ 4 ) 



The conclusion we can make is that there is a dramatic 
increase in the polarization cloud surrounding the exter- 
nal charges as one passes from the Higgs phase to the 
novel phase. From Fig. [13] and [14] we see that the coeffi- 
cient for the perimeter law for /i = .13 in the novel phase 
is approximately 9 times larger than that for [i = .2. In- 
deed, we can construct graphs analogous to Figs. [13] [14] 
for many values of and plot the parameter C, which 
is the slope of the line in the graph of — ln(W) versus 
L + T. We find a sharp cross over at the transition as 
shown in Fig. [15J indicating the almost ten-fold increase 
of the polarization cloud energy. 




2pi.q/e 



FIG. 14: (color online) Plot of -ln(W) for fx = .2. The 
dashed lines are a linear fit, we find — ln(VK) = 0, 001030 ± 
(3 x 10~ 7 )(L + T) + 8, 7 x 10~ 4 ± (5 x 10" 5 ) 



FIG. 16: (color online) Polyakov loop for various values of \i 
for a lattice 20xl00 2 . 



VI. SCALING STUDY 




FIG. 15: (color online) Plot of the parameter C as a function 
of fj, at fixed q. 



C. Polyakov loop 

At finite temperatures, one looks at the behaviour of 
the Polyakov loop operator, which is defined as the Wil- 
son loop variable taken along the entire (periodic) time 
direction N T for a fixed spatial position x. This is re- 
lated to the free energy of the system, F q , in the pres- 
ence of a single heavy quark by [23]: (P (x)) = e~^ Fq . In 
Fig. [l6j we see the behavior of the expectation value of 
the Polyakov loop operator (P) mirrors almost exactly 
the behaviour of the Wilson loop as in Fig. [6j The posi- 
tion of the transition has changed to a smaller value of 
and at a larger value of 2irq/eas should be expected, be- 
cause the temperature has been increased, corresponding 
to a Euclidean time direction of length 20. 



In any lattice simulation, it is important to use a lattice 
which is sufficiently large to eliminate finite size effects. 



In Fig. 17 the mean total length of loops as a function 
of \i is displayed for various lattice sizes; normalizing by 
the maximum possible length of loops Fig. |18| we see 
that it is independent of the lattice size, as expected. We 
find that the transition point occurs for /i rs 0.15 for all 
lattices larger than 10 3 , so the lattice size used in this 
study (100 3 ) is amply sufficient. 
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FIG. 17: (color online) Mean total length of loops {Lt) for 
various lattice sizes. 

In Fig. [l9j the density of finite loops versus the total 
density is illustrated, for various lattice sizes. The large 
error bars are only present for the smaller lattices sizes, 
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FIG. 18: (color online) Normalized mean total length of loops 
(Lt) I Lmax for various lattice sizes. 



FIG. 20: (color online) W(L, T) for fj, 
lattice sizes. 



0.13 for various 



10 3 , 20 3 , already at 40 3 we approach a consistent size in- 
dependent density as a function of /i. For lattices sizes 
greater than w 90 3 , the results are essentially indepen- 
dent of the lattice size. Finally in Fig.[20]the expectation 
value of two sizes of Wilson loop are displayed for various 
lattices sizes for /i = 0.13 as a function fo 2irq/e. Again, 
the results are clearly independent of the lattice size. 



0,25 




FIG. 19: (color online) Density of finite loops as a function 
of the total density for various lattice sizes. 

Therefore we conclude that using a lattice considerably 
larger than 40 3 , and especially the size 100 3 that was 
used for most of the simulations are perfectly adequate 
to remove all finite size effects. 



VII. DISCUSSION AND CONCLUSIONS 

Our results show numerical evidence for a novel phase 
in the phase diagram of the 3-d abelian Higgs model at 
the asymptotic boundary corresponding to strong (infi- 
nite) coupling in the spontaneously broken Higgs phase. 
At strong coupling the perturbative massive excitations, 
corresponding to the gauge boson and the neutral scalar, 
are completely decoupled. The only remaining particle 
is the vortex, which is an adjustable parameter. For a 
large mass of the vortex, the vacuum configuration is sat- 
urated by short loops of virtual vortex anti-vortex pairs. 
As this mass is lowered, the vacuum is filled with longer 
and longer virtual vortex anti-vortex pair loops. Finally 
at a critical value, there is a transition to a novel phase, 
in which the vortex loops reorganize into one effectively 
infinite loop in addition to a bath of smaller finite loops. 

In the Higgs phase, external charges are screened due 
to a polarization cloud which leads to a perimeter law 
for the Wilson loop, external charges are not confined. 
Smaller than a critical value for the mass of the vortices, 
the polarization cloud increases dramatically causing the 
energy shift as defined by the Wilson loop to increase 
9-10 fold. The Wilson loop behaviour however, remains 
the perimeter law, contrary to the behaviour that was 
surmised in [TU]. Since we have decoupled all perturba- 
tive excitations of the scalar field, including specifically 
charged excitations, it is in principle possible for the Wil- 
son loop to exhibit linear confinement. We explicitly find 
no dependence on the area for the Wilson loop. We find 
that the Wilson and Polyakov loop order parameters both 
vanish after a large enough value of the external charge, 
however this simply means that our lattice calculation is 
not able to resolve the details of its behaviour. 

The novel phase is characterized by the appearance of 
an effectively infinite vortex loop. The individual finite 
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vortex loops suddenly reorganize at the transition into 
one infinite loop and a gas of remaining finite loops, as a 
function of /i. The total length of the loops increases as a 
function of decreasing vortex mass, primarily through ap- 
pending to the infinite vortex loop. This novel phase was 
predicted by [10] and also in [IT]. In [11] the transition 
is described as the passage between the phases labelled 
VII and I. 

In [23] it was proven that there is no transition between 
the Higgs phase and the Coulomb phase, however, they 
look strictly at the compact version of the model. Wc 
are considering the non-compact model here, the analy- 
sis for non-existence of confinement does not apply. Un- 
fortunately, we do not find any evidence of confinement, 
even for fractionally charged external charges. Numerical 
work similar to ours was also undertaken by [25] , however 
they did not see our phase. There are two differences be- 
tween our simulation and theirs. First, we simulate only 
the effective model; hence, our results are valid only on 
the asymptotic boundary of the full phase diagram of the 
Abelian Higgs model whereas they [25] , simulate directly 
the gauge fields and scalar fields. Second, the limits taken 
are slightly different: we take A,e — > oo keeping A/e 2 
fixed, whereas they take first A — > oo and only afterwards 
do they take e — > oo. Hence we go to the corner of the 
phase diagram (in A and e space) along a particular fixed 
line, while they go up to the rectangular edge first (at 



A = oo) and then move along the edge (taking e — > oo) 
to the corner. It is clearly possible the two limiting pro- 
cedures do not commute. Furthermore, their numerical 
work was done on a 16 3 lattice, which, according to our 
scaling analysis, most probably will exhibit finite size ef- 
fects. The existence of the novel phase is expected to 
have important ramifications for the phase structure of 
the model in the presence of the Chern-Simons term |26) . 



VIII. ACKNOWLEDGMENTS 

We thank NSERC of Canada and the Center for 
Quantum Spacetime of Sogang University with grant 
number Rl 1-2005-021 for financial support. We thank 
Yvan Saint-Aubin, Jeong-Hyuck Park, and David Rid- 
out for useful discussions and we thank Jacques Richer, 
of the Reseau quebecois de calcul de haute performance 
(RQCHP), Montreal, for paramount help with the pro- 
gramming and its implementation. We also thank the 
RQCHP for allocating computer time. Finally we thank 
the Perimeter Institute for Theoretical Physics and the 
(Kavli) Institute of Theoretical Physics of the Chinese 
Academy of Sciences for hospitality while the paper was 
being revised. 



[1] P. W. Anderson, Phys. Rev. 130, 439 (1963). 

[2] P.C. Hohenberg, Phys. Rev 158, 383 (1967);N.D. Mer- 

min, H. Wagner, Phys. Rev. Lett. 17, 11331136 (1966); 

Sidney Coleman, Comm Math. Phys. 31, 259 (1973). 
[3] R. Dashen, B. Hasslacher and A. Neveu, Phys. Rev. D 

10, 4130 - 4138 (1974). 
[4] P. W. Higgs, Phys. Lett. 12, 132 (1964), ibid. Phys. Rev. 

Lett. 13, 508 (1964), ibid. Phys. Rev. 145, 1156 (1966); 

F. Englert and R. Brout, Phys. Rev. Lett. 13, 321 (1964). 
[5] S. Coleman and E. Weinberg, Phys. Rev. D7, 1888 

(1973). 

[6] A. Rajantie, Physica B255, 108 (1998); M. Chavel, Phys. 

Rev. Lett. B378 227 (1996); M. G. Amaral and M.E. Pol, 

Z. Phys. C44, . 
[7] A.M.Polyakov, Nucl. Phys. B120, 429 (1977). 
[8] H.B. Nielsen and P.Olesen, Nucl. Phys. B61, 45 (1973). 
[9] K. Wilson, Phys. Rev. D10, 2445 (1974). 
[10] S. Samuel, Nucl. Phys. B154, 62 (1979). 
[11] M. Einhorn and R. Savit, Phys. Rev. D19, 1198 (1979), 

see also T. Banks, R. Myerson and J. Kogut, Nucl. Phys. 

B129, 493, (1977). 
[12] M. Creutz, L. Jacobs, and C. Rebbi, Phys. Rev. Lett. 42, 

1390 (1979); Phys. Rev. D20, 1915 (1979). 
[13] Bollobas, Bela; Riordan, Oliver (2006), Percolation, 

Cambridge University Press, Grimmett, Geoffrey (1999), 

Percolation, Springer Verlag, 
[14] L. Jacobs and C. Rebbi, Phys. Rev. B , 19,4486,1979. 
[15] M. Sakellariadou, Nucl. Phys. B468, 319 (1996); M. 

Sakellariadou and A. Vilenkin, Phys. Rev. D37, 885 

(1988); A. G. Smith and A. Vilenkin, Phys. Rev. D36, 

990 (1987); T. Vachaspati and A. Vilenkin, Phys. Rev. 



D30, 2036 (1984). 

[16] R. Hagedorn, Nuovo Cim. A56, 1027-1057,1968. 

[17] J. J. Binney, N. J. Dowrick, A. J. Fisher and M. E. J. 
Newman, 'The theory of critical phenomena- an intro- 
duction to the renormalization group', Clarandon Press, 
Oxford, (1992). 

[18] D. P. Landau and K. Binder, 'A guide to Monte Carlo 
simulations in statistical physics', Cambridge University 
Press, (2000). 

[19] Madras, N. and Slade, G. 'The Self- Avoiding Walk', 
Birkhuser, Boston, (1993). 

[20] D. Mitchell and N. Turok, Phys. Rev. Lett. 58, 1577 
(1987), S. Frautschi, Phys. Rev.' D 3, 2821, (1971), R.D. 
Carlitz, Phys. Rev. D 5, 3231, (1972). 

[21] S. Coleman, 'Aspects of symmetry', Cambridge Univer- 
sity Press, (1999). 

[22] Bochner, Salomon, 'Harmonic analysis and the theory of 
probability', University of California Press, (1955). 

[23] H. J. Rothe, 'Lattice gauge theories - An introduction', 
World Scientific Publishing, (1992). 

[24] E. Fradkin and S. Shenker, Phys. Rev. D 19, 3682 (1979). 

[25] M. N. Chernodub, F. V. Gubar ev and M. I. Polikarpov, 
Phys. Lett. B 416, 379 (1998), [arXiv:hep-lat /970402l| , 
see also E. T. Akhmedov, M. N. Chernodub, M. I. Po- 
likarpov, and M. A. Zubkov, Phys. Rev. 53, 2087, (1996) 
and M. I. Polikarpov, U. J. Wiese and M. A. Zubkov, 
Phys. Lett. B309, 133, (1993). 

[26] P. W. Irw in and M. B. Paranjape, |arXiv:hep- 
th/9906021] . 



